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ABSTRACT 

Morphological and spectroscopic evidence suggest that shocks may affect the spatial 
and velocity distributions of gas in the narrow line region (NLR) and extended NLR of 
some active galaxies. It thus seemed plausible that shocks may also energize the NLR. 
The observed emission line ratios strongly favor photoionization as the heating source 
for this gas, but it is not clear whether the ionizing radiation is generated in the NLR 
by "photoionizing shocks" or whether the ionizing radiation originates at the central 
continuum source. Here I point out that shocks are highly inefficient in producing line 
emission. Shocks in the NLR can convert at most ~ 10~^ of the rest mass to ionizing 
radiation, compared with a maximum conversion efficiency of ~ 10~^ for the central 
continuum source. The required mass flow rate through "photoionizing shocks" in the 
NLR is thus a few orders of magnitude higher than the mass accretion rate required to 
power the NLR by the central continuum source. Since gravity appears to dominate the 
NLR cloud dynamics, shocks must lead to an inflow, and the implied high inflow rates 
can be ruled out in most active galaxies. 

NLR dynamics driven by a thermal wind or by some jet configurations may produce 
the mass flux through shocks required for photoionizing shocks to be viable, but the 
mass flux inward from the NLR must be kept ~ 100-1000 times smaller. 

Photoionizing shocks are a viable mechanism in very low luminosity active galaxies 
if they are highly sub-Eddington 10~^) and if they convert mass to radiation with a 
very low efficiency 10~^). 



Subject headings: galaxies: active-galaxies: Seyfert- shock waves 
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1. INTRODUCTION 



2. QUALITATIVE ESTIMATES 



Active Galactic Nuclei (AGNs) display emission 
lines originating from clouds in the narrow line re- 
gion (NLR) and the broad lines region (BLR; e.g. 
Netzcr 1990). The line emission could be powered ei- 
ther locally through cloud-cloud collisions (Daltabuit, 
Mac Alpine & Cox, 1978), or by absorption of ioniz- 
ing radiation from the continuum source at the cen- 
ter (e.g. Davidson & Netzer 1979). The time lagged 
response of the broad emission lines to UV contin- 
uum variability (e.g. Peterson 1993) provides nearly 
conclusive evidence that the BLR clouds are powered 
by the central continuum source. The same relation 
cannot be established for the NLR since reverbera- 
tion timescales are far too long. The observed line 
ratios in the NLR argue strongly in favor of pure pho- 
toionization (e.g. Ferland & Osterbrock 1986), which 
suggests that the NLR is also powered by the cen- 
tral continuum source. However, some models suggest 
that shocks due to cloud-cloud collisions may provide 
significant heating for the NLR, in addition to the 
central continuum source (e.g. Contini & Aldrovandi 
1983, Contini 1997). 

Renewed interest in shocks as a viable energy 
source came following the detailed calculations of Do- 
pita & Sutherland (1995, 1996) who have shown that 
fast shocks produce strong ionizing radiation which 
can be reprocessed into line emission. Thus, although 
the emission line spectrum from the NLR strongly 
suggests that clouds are powered by photoionization, 
rather than coUisional ionization, this power could be 
produced by shocks in the NLR, rather than by the 
central continuum source. Morphological and spec- 
troscopic evidence suggest that jet-cloud interactions 
may produce such shocks in the extended NLR of 
nearby AGNs (e.g. Cecil, Morse, & Veilleux, 1995; 
Capetti, Axon, & Macchetto 1997; Falcke, Wilson & 
Simpson 1998). For a comprehensive review and a 
critical assessment of the available observational ev- 
idence for photoionizing shocks in AGNs see Morse, 
Raymond & Wilson (1996) and Wilson (1997). 

The purpose of this note is to point out that 
shocks are highly inefficient, compared with the cen- 
tral AGN, in converting mass to ionizing radiation, 
and can thus be ruled out as a significant source of 
energy for the NLR for most AGNs. 



Shocks convert part of the directed kinetic energy 
of the colliding gas clouds to thermal energy. If all the 
kinetic energy of cloud-cloud collision in the NLR is 
converted to heat, all the heat is converted to ionizing 
radiation, and all the ionizing radiation is absorbed 
and converted to line emission, then the maximum 
fraction of rest mass converted to line emission would 
be Cshock = . The typical velocity of gas in 

the NLR is ~ 250 - 500 km s-i(sec §3), and thus 
the maximum efhciency of rest mass to line emission 
conversion is Cshock ~ 10^^. 

If the gas reaches the center of the active nucleus, 
about 10% of its rest mass is converted to heat, which 
may be radiated away and converted to line emission, 
yielding a maximum efficiency of eagn ~ 0.1. Thus, 
basic energy conservation implies that the photoion- 
izing shocks scenario requires a mass flow through 
NLR shocks which is a few orders of magnitude larger 
than required into the central photoionizing source. A 
more accurate estimate of the required accretion rates 
in both scenarios is given below. 

3. QUANTITATIVE ESTIMATES 

The [O III] A5007 hnc is typically the strongest 
narrow line in bright active galaxies (e.g. Ferland & 
Osterbrock 1986; Veilleux & Osterbrock 1987), and I 
therefore use it below as a measure of the NLR lumi- 
nosity. 

Dopita & Sutherland (1995) derived the following 
scaling relation for the H/? flux per unit area gener- 
ated in the cooling gas behind a fast photoionizing 
shock 

= 7.44 X lO^'^Wioo^n ergs cm"^ s"^ , 

where v = lOOtiioo km s~^ is the shock velocity and 
n is the ambient gas density. Additional line emission 

is produced in the unperturbed gas in front of the 
shock, the " precursor H II region" , which is given by 

lYid = 9.85 X lO^^iof " ergs cm'^ s"^ . 

The highest [O III] A5007/H/3 flux ratio found by Do- 
pita & Sutherland is ~ 13, and we therefore assume 
^[O III] ^ 13/h/3- The mass flux through the shocks is 

Mshock = A^nv, 

where Ag = i[o iii]/^[o iii] is the shock area, and 
-^[O III] is the intrinsic [O III] A5007 luminosity. Thus, 
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using the shock + precursor Hne emission, one gets 

M,hock > 2730i42t^roo^'(l + l-32t^roo'^)"' ^0 Y^"' , 

where L[o iii] = 10^^L42 erg s~^. 

The relation between the mass accreted by the 
AGN, Magn, and L^q uj] resulting from photoion- 
ization by the AGN is not unique. This relation de- 
pends on the accretion efficiency, the covering fac- 
tor of the NLR gas, the ionization parameter and 
density of the NLR gas, and the presence of dust 
within the NLR gas. Rather than deriving a the- 
oretical relation between Magn and L^q m], which 
would require knowledge of the parameters mentioned 
above, I use the observed L[o iii] versus Magn re- 
lation. The continuum i/L,^ at 5000 A is given 
by X[o ni] x 5000/£^W^[o ni], where EW^o iii] is the 
equivalent width of [O III] in units of A. The bolo- 
mctric luminosity is LboI ^ 10 x i^_Li^(5000 A) (e.g. 
fig. 7 in Laor & Draine 1993), and the accretion rate 
is Magn = -^Boi/eagnc^- These relations give 

AfAGN = 0.293U2EWioho\ Mq yr^^ , 

where EW[o in] = 'iOEW^o A, and eAGN = O.leo.i- 

The ratio Mghock/MAGN is a function of wioO) EW30 
and eo.i- Below we estimate typical values and dis- 
persions for these parameters. 

The [O III] FWHM data for 71 broad line active 
galaxies presented in Fig. 6a of Whittle (1985) gives 
(vwo) = 3.5^}'"', where the quoted range includes 2/3 
of the objects (equivalent to ±lcr for a normal distri- 
bution). Data on the EW[o iii] of 36 Seyfert 1 galax- 
ies in Ostcrbrock (1977) gives {EW[o iii]) = 39 ± 22, 
and the complete sample of Boroson & Green (1992) 
gives {EW[o III]) = 37^23 for 18 Seyfert 1 galaxies 
and {EW[o iii]) = '20tll for 69 quasars. We therefore 
adopt vioo ~ 3.5, and EW30 ~ 1 as typical values, 
which gives 



with a typical range of about a factor of two above 
and below. 

The value of eo.i can bo cstimatcid by comparing 
the time integrated density of quasar light with the 
current density of massive black holes, as determined 
by the black hole to bulge mass correlation suggested 
by Magorrian et al. (1997). This comparison provides 
a rather strict lower limit on AGNs time average ac- 
cretion efBciency of eo.i — 0.09 (see equations 1 & 2 in 



Hachnclt, Natarajan, & Rccs 1997 with h = 0.7 and 
Jb = 0.1). Larger values for the mean observed eo.i 
are obtained if not all bulges contain massive black 
holes, or if some of the black holes growth occurs in 
an unobserved phase with cq.i <C 1 (see discussion 
in Haehnelt et al.). Typical accretion scenarios give 
eo.i ^ 1- Thus, Mghock is at least about two, and 
more likely nearly three orders of magnitude larger 
than Magn- 

4. DISCUSSION 

4.1. Gravity Powered Shocks 

If the gas velocity in the NLR is dominated by 
gravity, as strongly suggested by observations (Whit- 
tle 1992a, 1992b), then the loss of kinetic energy in 
shocks necessarily implies a loss of angular momen- 
tum, which must lead to an inflow. If there is no 
accumulation of mass somewhere between the cen- 
tral black hole and the NLR then Afghock eventually 
becomes Magn- Having Mghock — -Magn requires 
eo.i ~ 10~^. However, the estimates of Haehnelt 
et al. mentioned above provide a lower limit of 
eo.i ~ 0.1 for the average efficiency. Thus, if the NLR 
of most AGNs is powered by shocks, then the mass 
flow through the NLR cannot accumulate in the cen- 
tral black hole. Can Mghock accumulate somewhere 
between the central black hole and the NLR? 

The black hole mass estimates of Magorrian et al. 
and of earlier studies of the cores of normal galax- 
ies (e.g. Kormendy & Richstone 1995) are based on 
stellar velocity dispersions in the central 10-100 pc. 
Thus, the limits on the accretion efficiency apply on 
the scale of the NLR as well. Stellar velocity disper- 
sions in a few nearby AGNs also indicate that the 
mass interior to the NLR is significantly smaller than 
expected from shock excitation (Oliva 1997). Thus, 
as pointed out by Oliva (1997) for the coronal line re- 
gion in NGC 1068, there is no place to hide the large 
-Mghock implied by the photoionizing shocks scenario 
for an astrophysically interesting period of time. 

Cloud-cloud collisions as a source of photoioniz- 
ing shocks are a viable mechanism only in low lu- 
minosity AGNs where L/Le <^ 10~^, such as the 
Galactic center (Narayan, Yi, & Mahadevan 1995), 
M87 (Reynolds et al. 1996, see a photoionizing shock 
model in Dopita et al. 1997), and possibly LINERS 
(Ho, Filippenko, & Sargent 1993), since such objects 
can maintain eo.i ^ 10~^, without accumulating ex- 
cessive mass, for a non negligible fraction of their 
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lifetime. The observed emission spectrum of LIN- 
ERS, however, appears to argue against photoionizing 
shocks (Ferland & Netzer 1983; Filippenko 1985; Ho, 
Filippenko, & Sargent 1993, 1996; Maoz et al. 1995). 

4.2. Wind Powered Shocks 

Shocks in the NLR would not necessarily lead to 
an inflow if the gas motion is not dominated by 
gravity. In particular, the shocks may be produced 
by interaction of NLR clouds with an outflowing 
wind. The clouds cannot be in pure radial outflow 
since that would result in mass ejection at a rate 
Afshock and mass and energy conservation requires a 
time average Mghock < -^agn, which can be main- 
tained only when eo.i ^ 10""^. One may imagine in- 
stead a "turbulent" velocity fleld where the required 
-^shock ^ 22OL42M0 yr~-^ (for uioo = 3.5) flows in 
and out within the NLR, with only 0.1-1% of this 
accretion rate being able to flow further inward to 
provide Magn- 

What can power such a wind? One possibility may 
be a radiation-pressure driven wind (e.g. on dust 
grains). The maximum momentum flux available in 
this case is Prad = Lbo\/c, or P^ad = Magncagnc 
The momentum flux in shocked gas is Pshock = 
AlshockV, and the maximum ratio between available 
and required momentum fluxes, 

Aad/Pshock = 0.1eo.ic/750eo.iw, 
is about an order of magnitude too small. 

Another possibility is a thermal-pressure driven 
wind, in which case momcntun is not conserved lo- 
cally, but only globaly. What can heat this wind? The 
wind velocity should be of the order of 250-500 km s~^ 
which corresponds to T ^ 1 — 3 x 10^ K. This temper- 
ature range is obtained naturally for Compton-heated 
winds in AGNs (Begelman, McKee, & Shields 1983; 
Mathews & Ferland 1987). The photoionizing energy 
flux in the NLR is E'phnphC, where i^ph ~ 30 eV is the 
mean photon energy, and riph is the photon density. 
The mechanical energy flux provided by the wind is 
Ep2nc, where iJp ~ 1 keV is the mean particle energy. 
The photoionizing/mechanical flux ratio is thus 13U, 
where U = riph/n is the wind ionization parameter. 
The Compton temperature is obtained for U > 100, 
and thus mechanical (i.e. shock) heating would be a 
negligible eff'ect compared with direct photoionization 
heating. 

Alternatively, a thermal wind may be mechani- 
cally heated by the dissipation of a jet kinetic en- 



ergy, as seen for example; in jet cocoons. The rate 
of energy deposition by the wind into the shocks is 
roughly ^Mshock?^^ = 3.2 x 10^^^^^ erg s"^ (for 
'-"loo = 3.5), which should be provided by the jet. 
Since LboI = 1-7 x 10'^^L42EW^q^ , the jet requires 
only ~ 0.2% of LboI- The energy deposited in the 
shocked gas should alos be observable as an extended 
non- variable source of continuum emission in the opti- 
cal to soft X-ray range. Additional spatially extended 
free- free emission atP^l — SxlO^K should come 
from the thermal wind, but the luminosity of this 
component could be well below 10*^ erg s~^. 

4.3. Jet Powered Shocks 

Although gravity appears to play the key role in 
the NLR, Whittle (1992a, 1992b) flnds evidence for 
excess velocity dispersions in a small fraction of AGNs 
(those with a high radio luminosity and linear radio 
morphology), which can be attributed to a jet-cloud 
interaction (see also Bickncll et al. 1997). Thus one 
can imagine a picture where clouds are swept up by 
the jet, and later fall back in, with nearly zero net 
accretion rate through the NLR, just as in the wind 
powered shocks picture. 

The jet starts at the center with velocity t;jet and 
mass flux Mjet, gradually entraining local gas until its 
mass flux increases to Mghock and its velocity drops 
to V. Mass, momentum, and energy conservation are 
used below to constrain the properties of such a jet. 
The jet motion is assumed to be balistic and any pos- 
sible interaction with a confining medium is neglected. 

Momentum flux conservation gives 

Mjet = f Mshock/7jot/3jetC , 

where /3jet = ?;jet/c and 7jet = (l-/3fgt)"^/^ (note that 
/^shock ^ 0.0017, 7shock ~ 1). Mass flux conservation 
requires Mjet < Magn which together with the above 
expression gives 

JjetPjetc/v > Mghock/MAGN • 

Using the Mghock/MACN relation obtained in §3 (note 
that this relation is independent of the presence of a 
jet), and v = 350 km s~^ gives "fjetl^jet > 0.875eo.i or 

/3jet > 0.875eo.i/y^l-h0.77e2 i . 

This is the minimum jet velocity which can provide 
the required momentum flux without exceeding the 
mass flux available through accretion. 
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The flux of kinetic energy in the jet must be smaller 
than the energy flux generated by the accretion, i.e. 

(7jet - l)MjetC^ < O.leo.lMAGNC^ • 

Substituting for Mjet from the momentum flux conser- 
vation, and using the Mghock/AfAGN relation obtained 
in §3 gives 

0.1eo.i7jet/3jetc/(7jet - 1)^ > 750eo.i , 

which for v = 350 km s~^ gives 

/?jet < 0.26 . 

This is the maximum jet velocity which can provide 
the required momentum flux without exceeding the 
energy flux available through accretion. The upper 
and lower limits on /3jet obtained above allow a solu- 
tion for /3jet only if eo.i < 0.26. Thus, for jet entrain- 
ment to be viable the jet must be subrelativistic, and 
the accretion efflciency must be fairly low. 

As the jet entrains ambient gas it slows down, 
while conserving momentum and losing kinetic en- 
ergy. The lost kinetic energy is converted to internal 
energy of the entrained gas and some of it will be 
radiated away. For example, in a jet with a solid 
angle of ^l/An = 0.01, reaching a final vioo = 3.5, 
i.e. Mghock = 22OL42MQ yr~^, at r = 10 pc, and 
T = lO^'K, the free-free cooling time is ^ 6 yr, which 
is well below the dynamical time scale, r/v ~ 3000 yr. 
The jet cooling will produce an extended continuum 
source, in addition to that within the shocked gas it- 
self (§4.2). This continuum can be quite intense. The 
initial jet kinetic energy is L™^ ~ ^||Z/bo1) the ac- 
cretion efficiency is eo.i > 0-1 (§4-1), which gives a 
minimum jet velocity of /3jet — 0.087 (from momen- 
tum conservation). Since the final jet kinetic energy 
is only L^^^^ ~ 0.005Lbo1, the dissipated kinetic en- 
ergy in the entrained gas is large ~ 0.3Lbo1- A sig- 
nificant fraction of this internal energy is likely to 
be radiated through thermal and nonthermal contin- 
uum emission somewhere in the X-ray to radio bands. 
Thus, observational constraints on the amount of spa- 
tially extended continuum emission in nearby AGNs 
can be used to constrain the feasibility of the jet-cloud 
interaction scenario. Wilson, Ward & Haniff (1987) 
find that the extended radio power in nearby AGNs is 
< 10~^Lnlr- Comparable constraints in other bands 
would rule out jet entrained gas as the source of power 
for the NLR in most AGNs. 



If the jet is still confined as it entrains ambient 
gas, if cooling is negligible, and if the jet cross section 
increases outward, then the random particle motion 
can be converted back to directed motion. The j(;t 
basically converts kinetic energy to internal energy 
and back to kinetic energy, conserving kinetic energy 
and not conserving momentum (due to the confining 
medium). The limits obtained above will not be valid 
(see Phinney 1983 for a comprehensive discussion), 
and the jet just needs to provide the kinetic energy 
deposited in the shocks (0.002Lboi), as in the case of 
the thermal pressure driven wind discussed above. 

A jet may also power the NLR by directly dumping 
its kinetic energy into gas clouds in the NLR, rather 
than by inducing gas motion which results in pho- 
toionizing shocks, as assumed above. However, in 
this case the line emission would result from thermal 
excitation, rather than photoionization, which is not 
consistent with the observed spectrum from the NLR. 

5. CONCLUSIONS 

Photoionizing shocks due to cloud-cloud collisions 
in a gravity dominated NLR cannot provide the 
source of energy for narrow line emission in most 
AGNs. 

NLR dynamics dominated by forces other than 
gravity may allow the high mass flux through shocks 
required for photoionizing shocks to be viable. The 
clouds velocities can be produced by a wind-cloud in- 
teraction, which requires a mechanically heated ther- 
mal wind. Alternatively, the clouds velocities may 
be dominated by a jet-cloud interaction. This either 
implies a strong extended continuum source, which 
is not yet observed, or otherwise specific jet config- 
urations. In any non-gravitational cloud dynamics 
picture the mass flux inward of the NLR must be 
kept ~ 1 00-1000 times smaller than through the NLR 
shocks. 

Photoionizing shocks could be relevant in very low 
luminosity AGNs, such as LINERs, if these objects 
are both highly sub Eddington {L/Le <, lO^**) and 
if they convert mass to radiation with a very low ef- 
flciency (eAGN <, 10-^). 
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